PUPT-2364 



Conformal symmetry and the 
Balitsky-Kovchegov equation 



Steven S. Gubser 
Joseph Henry Laboratories, Princeton University, Princeton, NJ 08544 



Abstract 

Solutions to the Balitsky-Kovchegov equation are considered which respect an 5*0(3) 
subgroup of the conformal group. The symmetry dictates a specific dependence of the 
saturation scale on the impact parameter. Applications to deep inelastic scattering are 
considered. 
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1 Introduction 

The Balitsky-Kovchegov (BK) equation [H El E] is a central tool for understanding the initial 
conditions in hadronic collisions in situations where transverse gluon density approaches 
unitarity limits. At leading order in large N c and small a s , the equation reads 
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dS(r 1 ,r 2 ; Y) _ 
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a s = \ 2 ) 

7T 

is essentially the 't Hooft coupling, and r±, r 2 , and z are positions in the transverse plane, 
conveniently thought of as complex numbers. The function 

S(T U rxY) = ^(X,Trt(n)U(r 2 )) Y (3) 

is the correlator of fundamental Wilson lines in the presence of a hadronic target, evaluated 
at rapidity Y relative to the target. The range of S is [0, 1], and for fixed Y one usually 
requires S , (r 1 ,r 2 ;F) — > 1 as r\ — > r 2 , while S(r 1 ,r 2 ;Y) — > as r x moves far away from r 2 . 
The quantity 

M(n,r 2 ;Y) = l-S( ri ,r 2 ;Y) (4) 

is the amplitude for a quark-antiquark dipole, characterized by position t\ for the quark and 
r 2 for the anti-quark, to scatter off a hadronic target in a collision where the difference in 
rapidity between the dipole and the target is Y . A common application of the BK equation 
(though far from the only one) is to study deep inelastic scattering (DIS), in particular the 
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total cross-section of a virtual photon off of a proton. In this application, the scattering 
is presumed to be dominated by the overlap of the photon with a quark-antiquark dipole, 
whose likelihood to interact with the hadronic target is then described by J\f(ri,r2',Y). 

A standard simplification is to require that S depends on r\ and r 2 only through \ri — r 2 |. 
Physically, this means that the dipole size is assumed to be much smaller than the typical 
size of the distribution of hadronic matter with which it collides. This paper is focused on 
the observation that an equally valid simplification of the leading order equation is to 
require that 5 depends on r x and r 2 only through the distance function 

dq ^ r2) = m Z 21 'mn^ 21 m ' ( 5 ) 

a/(1 + g 2 |ri| 2 )(l + q 2 \r 2 \ 2 ) 

where q is an arbitrary parameter with dimensions of inverse length. 1/q is essentially the 
size of the hadronic target. The reason that solutions to of the form 

S(r 1 ,r 2 ;Y) = S q (d q (r 1 ,r 2 );Y) (6) 

must exist is that the functional form on the right hand side is the one that respects an 
5*0(3) subgroup of the symmetry group P5L(2, C) of the measure factor d 2 z \ ri _l^2-z\ 2 • 
This is essentially the same 5*0(3) subgroup as the one used in [4, 5J to find generalizations 
of Bjorken flow with finite transverse extent. Because two out of the three generators of 
50(3) are conformal Killing vectors rather than true Killing vectors, solutions of the form 
(|6| will be modified by the breaking of conformal invariance, which occurs at next-to-leading 
order (NLO) in the BK equation through the running of a s . The ansatz is similar in 
spirit to the work of [U] on solutions to the Banfi-Marchesini-Smye equation preserving a 
50(2, 1) subgroup of the conformal group. 

The organization of the rest of this paper is as follows. In section [2] I briefly review 
some of the standard lore on solutions of the leading-order BK equation, including the 
key concepts of saturation scale and geometric scaling. In section [3] I explain the 50(3) 
symmetry more fully and show how this symmetry dictates the functional form of the distance 
function d q (x,y). I will also demonstrate that the 50(3) symmetry, together with some 
weak additional assumptions, dictates the following dependence of the saturation scale on 
the impact parameter b: 

Qmax(Y) 

provided Qf &yi {Y) ^> q, which occurs at large enough rapidity. In section |4| I explain how 
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S(r%, r 2 ; Y) is converted into a prediction of DIS cross-sections, and point out that although 
the form ([7| is satisfactory at large Q, it leads to a total cross-section <7 7 * p that increases 
somewhat too quickly as Q becomes small. This mismatch is due to the failure of conformally 
invariant solutions to capture the rapid falloff of energy density at sufficiently large transverse 
radius. Conformal symmetry has also been used in [7] , along with some additional dynamical 
intuition from AdS/CFT, to study DIS cross-sections at small Bjorken x. 



2 Some well-known properties of the BK equation 

Reviews of the large literature on the BK equation and related concepts include [5J [91 [10] • 
Earlier literature, for example [TT1, IT2"]. shows that many aspects of the dynamics were under- 
stood before the BK equation was introduced in its modern form. The remarks in this section 
are intended to provide non-specialists with some perspective on how the BK equation is 
used. 

It is often convenient to replace r\ and r 2 by the combinations 

r = r\ — f2 and b = — ^ — - . (8) 

For consistency with the notation of Q, r and b are most properly regarded as complex 
numbers. However, when no ambiguity seems likely, I will abuse notation by allowing r and 
b to denote the magnitudes \r\ and |6|, respectively. In particular, throughout this section I 
will restrict attention to solutions to the BK equation which depend only on |r|, and I will 
use r to denote this positive real quantity. That is, 

S( ri ,r 2 ;Y) = S(r;Y). (9) 

An important and well-studied class of solutions is those which tend toward a scaling 
form at large Y: 

S(r;Y) -»■ S(r) as Y ->• oo , (10) 

where the scaling variable is 

f = Q s (Y)r (11) 

and Q S {Y) grows exponentially with Y: 

Q S (Y) = Q e^ aY where c w 4.88 . (12) 
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The estimate of c quoted here is from [T3J Q3] , where it is also argued that there are non- 



exponential corrections to (12). Q 8 (Y) is called the saturation scale. Intuitively, 1/Q S is 
the characteristic transverse size of gluons whose occupation numbers become comparable 
to unity. A typical value for Q s in the description of deep inelastic scattering from HERA 
[TBI \W\ \TT\ [T5] is 1 GeV [TJ] , though to assess precisely what this value means one must 
obviously have a specific form of S(f) in mind. Dependence of S(r; Y) and related quantities 
only on f (or mainly on f) is known as geometric scaling [19} 120] (see however the discussion 
following ( |38| ) for a refinement of what geometric scaling means when referring to DIS data). 
Intuitively, geometric scaling says that the scattering amplitude at rapidity Y for a dipole of 
transverse size r depends only on the number of gluons at the relevant rapidity in the part 



of the target sampled by the dipole, which is proportional to r 2 Q s (Y) 2 



It has been argued (21] that the approach (10) to geometric scaling is a universal feature of 



all solutions S(r;Y) to the BK equation satisfying certainly physically motivated conditions, 
in particular |1 — S(r;Y)\ <C r 2lc at some fixed Y for small r, where^j] 

7c ^ 0.628. (15) 



Qualitatively, the claim being made in ( 10 ) is that hadronic collisions have universal behavior 
at high energies, independent of the precise nature of the colliding particles. 

A certain amount is known about the scaling function S(r) that controls the asymptotic 
large Y behavior: 

• For r<l, one has S(r) = 1 + Cif 2lc logf + . . ., where C\ > is a constant [21J. 

• For f ~ 1, the numerical results of [22] fit well to the form S(f) = 1/(1 + f /f*) 2 , where 
r* is a constant. 

• For f > 1, one has S(f) = C*2e~c lo s 2r + . . ., where C2 > is a constant and c is the 



same constant as the one appearing in (12). 



^or the analysis of geometric scaling, it is convenient to start by defining 

AT(k;Y) = J d 2 re- ikr l ~^ Y) , (13) 
where k and r are regarded as vectors in R 2 . Then, remarkably, Q can be shown to be equivalent to 

^=a sX (-d L W-alAf\ (14) 



where L = logfc 2 and xil) — 2"0(1) — ^(7) — "0(1 — 7)- Here ip(z) = d\ogY (z) / dz . The exponent -f c in (15) 
is defined as the solution to the equation 7^(7) — x(l)- 
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However, it is far from clear that large Y asymptotics of the solutions of the leading-order BK 
equation are directly relevant for phenomenology. Parton fraction x ~ 10~ 4 , which approxi- 
mately characterizes the HERA deep inelastic scattering data as well as mid-rapidity LHC 
phenomena, corresponds to Y — — log a; = 9.2, whereas numerical investigations including 
[22] show considerable evolution of scaling forms even for Y > 20. So initial conditions at 
some modest value of Y are important. One favorite is the Golec-Biernat-Wiisthoff model 

m 

S GB w(r; Y ) = e~^> 2 , (16) 



where Y Q and Q s0 are parameters. 

An additional challenge is that with reasonable values of the coupling, for example d s 



0.2, the growth (12) is substantially more rapid than a comparison with DIS data supports. 
This situation is substantially improved by the inclusion of next-to-leading-order effects: see 
for example 



3 Exploiting conformal symmetry 

I want to consider distributions of hadronic matter which are non-uniform in the transverse 
plane. In the literature on the BK equation, this is relatively unexplored territory. In [8] 
one can find a form for the dipole scattering amplitude based on treating a nucleus as a 
sphere of uniform density. In [24J one can find calculations related to the Froissart bound 
based on the assumption of a factorized form N{r\,r%) = Af(r)S(b). Most relevant for 
our investigations, in [25] one can find the expression (|7|) as an ansatz for the dependence 
of the saturation scale on impact parameter. Earlier related work includes [2"fi"| ¥T7\ [2H], 
which consider how diffractive processes constrain the impact parameter dependence of dipole 
scattering amplitudes; and [29J, in which solutions to the BK equation with explicit impact 
parameter dependence were investigated numerically. 

In general, it's hard work to study the BK equation when S has some general dependence 
on complex r\ and r 2 rather than dependence only on r = \r\— r%\, simply because the general 
form S(rx,r2',Y) has five real independent variables, as compared to two independent real 
variables in S(r; Y). Dependence of S only on r = \r\ — r%\ and Y amounts to the imposition 
of ISO(2) symmetry. The action of ISO(2) on a complex number z (describing, as usual, a 
position in the transverse plane) is 

z^az + (3 (17) 
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where \a\ = 1 and is any complex number. Imposing ISO{2) symmetry on a function 
S(ri, r 2 ; Y) means that we demand that S(ri, r 2 ; Y) is invariant if we send ri — > ar\ + /3 and 
r 2 — > ar 2 + (3 at the same time (and leave Y unchanged). 

In jl] I showed in the context of relativistic hydrodynamics that an interesting solution- 
generating technique was to replace invariance under 750(2) by invariance under a group 
S0(3) of conformal isometries of R 3 ' . I will exploit the same strategy here. The reason it 
works is that the measure and kernel of the leading-order BK equation is invariant under 
the group PSL(2, C) of linear fractional transformations (LFTs)J^] which act on z by 

az + 8 , , 

jz + d 

where a, (3, 7, and S are complex. A simple way to check this invariance is to consider only 
the holomorphic half of the measure and kernel, namely the one-form 

77 — r 2 
U (77 - z)(r 2 - z) Z ' 

The claim is that u is invariant under all linear fractional transformations, where we un- 
derstand that the LFT is supposed to act on 77, r 2 , and z. To check this claim, first note 



that it is obvious for linear transformations (17) (even when a is not unimodular), and that 



the full group PSL(2, C) is generated by linear transformations together with the inver- 
sion z — > — 1/z. So the only thing to check is that u is invariant under z — > —1/z (with 
r% — > — 1/ri and r 2 — > — l/r 2 at the same time). This is a straightforward exercise. Con- 
formal invariance of the BFKL kernel has been thoroughly appreciated in the literature (see 
for example [12J ) , and the invariance of the BK equation that I have just described has been 
noted explicitly, for example in |29j . 

Just as J5'0(2)-invariance requires that S should depend on ri and r 2 only through 
the combination \r\ — r 2 |, so one would expect invariance under some other subgroup G G 
PSL(2, C) to require that S should depend on r\ and r 2 only through some other simple 
combination. If such a symmetry is imposed on S(r 1 ,r 2 ;Y) at some fixed Y, then the 
PSL(2, C) invariance of the leading-order BK equation ensures that it persists at all Y. 
The obvious example to start with is an SO (3) subgroup which includes rotations around 
the beam axis (that is, rotations of the complex plane preserving the origin). There is a 



2 SL(2, C) is the group of matrices where a, j3, 7, and 8 are complex with a5— /3j = 1. PSL(2, C) 

is SL(2, C) quotiented by the Z2 which simultaneously flips the sign of a, f3, 7, and S. LFTs of the form 



(18 1 are in unique correspondence with elements of PSL(2, C). 
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one-parameter family of such subgroups, characterized by a length scale which I will denote 
1/q, following p[]. Besides rotations preserving z = 0, SO(3) q includes the LFTs 



qz cos % + sin ^ 

^ -> -V 2 , , (20) 
— qz sin 2 + cos | 

which preserve the points 2; = ±i/q. The parameter G [0, 27r) uniquely labels these 
conformal isometries of the complex plane. The rest of SO(3) q can be obtained by composing 



rotations preserving z = with transformations of the form (20). It is useful to recall that 
there is a conformal map, the so-called stereographic projection, which maps the complex 
plane to the unit sphere S 2 . Let's parametrize the S 2 by polar angles (9, 0) such that 9 runs 
from at the north pole to tt at the south pole, and <ft runs from to 2tt and labels longitude. 
Then the stereographic map is 

Q 

\qz\=tsxi- axgz = (p, (21) 



and the transformations (20) can be understood as ordinary rotations of S 2 which preserve 
two opposite points located on the equator. 

The SO(3) q which I have described is essentially the same as the one in [?], and its 
original motivation came from collisions of gravitational shocks in AdS$ |30j [31] |^] More 
precisely, each element of SO(3) q is a conformal isometry of the complex plane which can be 
extended to a conformal isometry of R 3 ' 1 , and this extended isometry is an element of the 
group SO(3) q picked out in [1]. This is most easily seen at the level of derivative operators 



representing the generators of SO(3) q . The infinitesimal generator of transformations (20) 



is 



( ± = (l + q 2 z 2 )d z + (l + q 2 z 2 )d z . (22) 
If we write z = x 1 + ix 2 , then 



C ~ — + q 2 

dx 1 



(23) 



interestingly, there is a literature on comparing AdS/CFT calculations involving gravitational shocks 
to DIS data: see for example [35]. This work is, however, largely orthogonal to the present work, because 
the focus here is on the use of SO(3) q as an organizing principle for dependence of scattering amplitudes on 
impact parameter, whereas the focus in [32] is on the AdS/CFT prescription for computing Wilson loops. 
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Now recall the form of the symmetry generator used in [I] to define SO(3) q : 

d 



^ dx 1 ~^ ^ 



d d 
dxv ^ dx 1 



(24) 



where x^ = (t, x 1 , x 2 , x 3 ) and x M = (— t, x 1 , x 2 , x 3 ). Evidently, (± = ( provided t 2 = (x 3 ) 2 . 
In other words, (± on the light-front locus t 2 = (x 3 ) 2 extends to £ on all of R 3 ' 1 . 

Having established the main properties of the SO(3) q symmetry, let's ask next what real 
function of complex r\ and r 2 is invariant under it. That there should be essentially one such 
function makes sense from the point of view of counting real parameters: T\ and r 2 together 
comprise four real parameters, and each of the three generators of SO(3) q can be thought of 
as constraining one. From a more geometrical point of view, it's obvious that there should 
be essentially only one S , 0(3) g -invariant interval, because on the unit sphere we know what 
it is: geodesic distance between the two points in question. More convenient than geodesic 
distance for many purposes is chordal distance: if we embed S 2 into R 3 in the standard way 
(i.e. as the unit sphere centered on the origin), then the chordal distance between two points 



on S is the length of the line segment connecting them in R . Using the definition (21) of 
the stereographic map, one may straightforwardly check that the chordal distance between 
the stereographic images of two points r\ and r 2 is 2qd q (ri, r 2 ), where d q (ri,r 2 ) is the distance 
function defined in ([5]). It is possible to show directly that d q (ri,r 2 ) is unchanged when r\ 



and r 2 are replaced by S0(3) q images of themselves, for instance under the LFT map (20). 

The discussion so far makes clear that there are solutions of the BK equation that respect 
SO(3) q invariance, and that they take the form 

S(r 1 ,r 2 ;Y) = S q (d q (r 1 ,r 2 );Y). (25) 

It remains to calculate S q (d q ;Y). This is hard, in the sense that no analytical results are 
likely to be available. However, the physically interesting regime to consider is Q s ^> q: that 
is, the saturation scale is much smaller than the overall size of the target. Then there is a 
simple way of generating approximate solutions: If S(r; Y) is an ISO(2) symmetric solution 
which vanishes quickly away from small r = |ri — r 2 |, then one need only replace r — > d q (ri, r 2 ) 
in order to obtain an 5*0(3)-invariant function which is an approximate solution to Q. The 
reason this works is that the 5*0(3)-invariant function is significantly positive only when 
qd q (ri, r 2 ) 1. So it doesn't "notice" the curvature of S 2 , and acts the way solutions do on 
R 2 . As a special case, we can start from an J5'0(2)-invariant scaling solution S(r; Y) = S(f). 
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The corresponding approximate S'0(3)-invariant solution is 



S(n,r 2 ;Y) = S(Q™(Y)d q ( n ,r 2 )). (26) 
In the JS'0(2)-invariant solution, the quantity Q™ ax (Y) is exactly the Q S (Y) that came up 



in (11 ) and (12); but I have renamed it in (26) in order to make clear its new meaning in the 



5'0(3)-invariant solution as the maximum of the saturation scale over the transverse plane. 
When qd q (ri,r 2 ) <C 1, one can make the further approximation 

\ r l — r 2\ \ r \ 

aq{r 1 ,r 2 ) « - 2 = — — . (27) 

I | Q 2\ n+r-2 r 1 + g 2 2 v 7 



Plugging (27) into (26), one finds 

S'(r 1 ,r 2 ;y)« i 5(g s (&;y)r) (28) 

where 

O max (Y) 

is the saturation scale at a position z in the transverse plane. The form ([7]) has been 
considered previously [25] (see also |33j), but now we can see precisely what its distinguishing 
feature is: It is the S'0(3) g -invariant choice for the dependence of the saturation scale on 
transverse position. 



To derive (29) I used the assumption of geometric scaling. However, much weaker as- 



sumptions suffice. Assume we have SO(3)-invariance, so that (25) holds for some function 
S q (d q ; Y). Assume further that for fixed d q > 0, one has S q (d q ; Y) — > as Y — > oo. Lastly, 
assume that for any fixed value of Y, S q (d q ; Y) is a monotonically decreasing function of d q . 
Let's define Q™ ax (^) through the equation 

s '(orW r ) = K ' (30) 

where k G (0, 1) is some pre-specified number. Q™ ax (F) is unique because of the monotonicity 
property. Q™ ax (Y) exists for sufficiently large Y because S q — > as Y — > oo for d q > 0Q 
Q™ ax (F) — > oo as Y — > oo, again because S q — > as Y — > oo away from d q = 0. The regime 



4 This is a non-trivial observation because S q does not in general reach arbitrarily small values. If we 



specify a value of k that is not particularly small, say k — 1/2, then Qf ax (Y) as defined in (30) exists for 
all Y above a modest threshold. 
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we are interested in is Q™ ax {Y) 3> <?• in this regime, up to ambiguities suppressed by powers 
of q/Qf^iY), we can define the saturation scale Q s {b; Y) as the solution the equation 

d9 ( &_ ^' H 2^) = « 7 )- (31) 



Using the approximation (27) to simplify the left-hand side of (31), we arrive immediately 
at 0. 

4 Deep inelastic scattering 

I will first briefly review how the amplitude S(ri,r 2 ;Y) is used to calculate certain cross- 
sections in deep inelastic scattering. For a more complete account, see for example [T9l 18]. 
A key quantity for comparison with experiment is the total cross-section cr 7 » p of a virtual 
photon, denoted 7*, to scatter off the hadronic target, which at HERA is a proton. Standard 
calculations, summarized for example in [9J, lead to the formulas 

a rp (x, Q 2 ) = a T (x, Q 2 ) + a L (x 7 Q 2 ) (32) 

where 



a T , L (x,Q 2 ) = j d 2 n J d 2 r 2 J dv\y TtL (r,v,Q 2 )\ 2 Af( ri ,r 2 ;Y) , (33) 
and I have again used the convenient notation r = \ri — r 2 \. Also, 



l^l 2 = E4 + (1 - «fl Q)Kr{Q f rf + m 2 f K (Q f r) 2 ) 
71 f 

^ = 3^n £ &2 (4Q _ vfKo{Qfrf) 
71 f 



(34) 



where Kq and K\ are modified Bessel functions. Also, 

Q) = v{\ - v)Q 2 + m) . (35) 

\Pt and are the wave functions for splitting the transverse and longitudinal polarizations 
of the virtual photon into a quark at transverse position r\ and an anti-quark at transverse 
position r 2 , with r = \r± — r 2 \ as usual, v is the momentum fraction of the photon carried 
by the quark. The subscript / in the sums indicates the quark species, and rrif and Cf are 
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the mass and charge of each quark. The dependence of |^t,l| 2 on r\ and r 2 only through r 
amounts to the assumption that the photon is equally likely to be anywhere in the transverse 
plane. Q 2 is the photon's virtuality (positive for a spacelike photon), and Y = log - where 
x is Bjorken's variable. 

A well-known analysis [THl ED] adopts the following simplified form of (33): 

<xt,l(x,Q 2 ) = a J d 2 r J dv \^ L (r,v,Q 2 )\ 2 Af GBW (r;Y) , (36) 
where a is a constant with units of area and 

AWr;y) = l- e -^ (y)V . (37) 



In the limit m/ — > 0, the geometric scaling property of (37) translates into a dependence of 
&t,l{ x iQ 2 ) on Y and Q 2 only through the dimensionless ratio 

- - Q2 (38) 



Q.0O* ' 

Geometric scaling in DIS data refers, most properly, to the dependence of <jy* p (x, Q 2 ) and 
related quantities only on r. Evidently, this can only work if Q a (Y) is chosen properly, and 
a good fit to the data is achieved using 

Q S (Y) = (lGeV)e ai4 ( y - 8 - 1 ) (39) 

and o"o = 23 fib, with a sum over the three flavors u, d, s with mf — 140 MeV (a value 
small enough so that it only slightly modifies geometric scaling) [TH]. The simplification 



(36) leaves something to be desired, because instead of integrating a scattering amplitude 
jV(ri, T2] Y) over the positions of both the quark and the anti-quark, one starts by assuming 
translational invariance M{ri,T2]Y) = Af(r;Y). Translational invariance of the hadronic 



target automatically makes all cross-sections infinite, and this is avoided in (36) by replacing 
the integration over impact parameter by the factor (T - 



To understand the main implications of SO (3)-invariant scattering amplitudes (26) for 
DIS, let's consider the following ansatz: 

M{r l) r 2] Y) = l-e'\ Q ^ 2r \ (40) 
where r = \r\ — r%\ and b = I Oj±ra I _ j p[ an ^ Q use ^ e expression (TH) for Q s (b; Y), but to start 
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with let's consider a general dependence of Q s on the magnitude of b and on Y. A useful 
simplification is to note that when m/ = 0, 

f\ V (|» T (r,.,qT + |» l( r,.,^n „ ^ (Qr)2 ; Wr)4/4 . (4!) 
This form, though approximate, is convenient because it allows all but one integral in the 



expression for a 7 * p to be performed explicitly: starting from (32)-(34) and (40), one finds 



a 



em Q 2 / bdb I rdri———— = ^ I bdbG[-if—) (42) 



TP 



r °° i l-e-\ Q ^ 2r2 4a em r^^^f Q 2 



o Jo 



(Qr) 2 + (Qr) 4 /4 3 J \Q s (b;Y)* 



where 

G(9) = 7-log0-e 1/e ReEi(-l/^). (43) 
Here Ei is the exponential integral function, and 7 is Euler's constant. Let's define 

r==^— and u(b;Y) = Q l {b]Y) , (44) 
Q.On 2 ^ ^ Q S (Y) 1 J 



so that (42) may be re-expressed as 



ffrp ^^ G (_L_). (45) 



In (44), Q S (Y) is an average or characteristic value of Q s , whose precise definition we do not 
need to specify at this stage. In order for <t 7 * p to depend on Q 2 and Y only through the 
combination r defined in ( |44| ), u must be a function only of b, not of Y. For S'0(3)-invariant 
solutions, according to ([7]), we have 

1 + g 2 6 2 ' V ; 

so u is independent of Y precisely if Q s is a fixed fraction of Q™ ax - 

It is easy to see that G(8) is a monotonically decreasing function with 

- log 9 + . . . for 6 < 1 
G{9) = I , (47) 
-log0 + ... for6'>l. 

where . . . indicates subleading terms. As a rough estimate, in the expansions (47), we may 



12 



in replace log# by —1 for 9 <C 1 and +1 for 9 ^> 1, and split up the last integral in (42) as 



a rp (x,Q 2 ) 



bdb 



°° UJU Q s (h;Y) 2 



(4f 



where 6* is determined implicitly by the equation 



Q S (K;Y)=Q. 



(49) 



If Q is large, then (49) has no solutions, and one must effectively set b* = in (48). Then 



one finds <r 7 » p oc 1/r where r is defined as in (44), up to logarithmic corrections due to the 



coarseness of the approximations in (48). If Q is small, then 6* is large, and the first term 
in square brackets in (48) dominates. In this situation, ct 7 . p oc b 2 , again with logarithmic 



corrections. In summary, 



&ry* p OC 



(50) 



for t > 1 
for r < 1 . 

The analysis so far has been for a general dependence of Q s on impact parameter b, 



provided the ratio u defined in (44) depends only on b, not Y. Evidently, the large r 



behavior <r 7 * p oc 1/r is universal, and, pleasingly, this behavior is a good fit to the data 



[19. 20 1 . According to (49) and (50). the small r behavior of (J>y*p IS 8b probe of the behavior 
of Q s (b; Y) at large b — assuming, as seems reasonable, that Q s (b; Y) decreases monotonically 
to (or to values below A QCD ) as b increases. If we assume the form ffh for Q s (b; Y), then 
(49) can be solved to give <t 7 * p oc 1/y/r for small r. This is a faster increase with decreasing 



r than is supported by the data: Referring to Figure 2 of [20], we see that as r decreases, 
\/r(T 7 .p saturates to a maximum value near r = 1 and then decreases roughly like a small 
positive power of r to as far as the data reaches, roughly r = 2 x 10~ 3 (see also figure [l] 
below). This decrease is consistent with a dependence <t 7 * p oc 1/\/t, though the width of 
the available region of r is narrow enough that a firm conclusion about the functional form 
that describes it should not be drawn. 

We should not be too surprised that conformal invariance predicts a faster increase in 
o" 7 * p at small r than data supports: this increase is due to the relatively slow decrease of 
Q s (b;Y) at large b exhibited by ([7]). Confinement indicates that Q s should decrease more 
rapidly than 1/b 2 once one reaches a sufficiently large value of b. 



A weakness of the analysis I have presented is that the form (40) with Q s (b;Y) given 
by (J7]) is only an approximate presentation of S'0(3)-invariance, valid when r < 1/g. A 
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more precise expression of SO (3)-invariance would be to use (26); however this makes the 
integrations quite a bit more complicated. There is a way of seeing that this weakness is 
not a serious problem, as follows. The regime of parameters corresponding to HERA data 
is Q between 0.2 GeV and 20 GeV, with Q s on the order of 1 GeV. On the other hand, q 
should be chosen to be approximately the inverse radius of the proton, i.e. 0.2 GeV. Thus 
Q ^ q, and we can expect the main contribution to the integrals defining <r 7 * p to come from 
r 1/q. In short, in the physically interesting regime, we are seldom far from satisfying 
the desired inequality r <C 1/q, so the use of §f§ is sufficient in order to understand the 
qualitative features of the dependence of <r 7 . p on r. 



Another issue is that (40) does not solve the BK equation. Clearly, an improved analysis 
would be desirable. However, I expect that the essential feature of excessive growth of <t 7 » p 
at small r is a durable consequence of conformal symmetry, since it is related to the absence 
of confinement. 



Cutting off the integration over b in (42) at some maximum value 6 max approximately 
equal to the size of the proton improves the fit to measurements of the total cross-section 
cx 7 » p , including the decrease in a/™Vp as r decreases well below 1. This is not too surprising 
since the cutoff in b is a crude implementation of the effects of confinement. For comparison 
to data, I used 



dm /""max / 

a, 



TP 







bdbG ^ ] (51) 



with u(b) given as in (46) and G(6) derived as in (42) from an improved approximation to 



the integrated photon wave function, namely 

f\ mr,,,QV + \n(r,v,Q^) „ (l + j^j^) (52) 

where s = Qr. Calculations were performed with three massless flavors and a em = 1/133. 
It turns out that a good fit for be achieved over a significant range of choices for 

q with a cutoff 6 max = 0.69 fmj^] provided that one allows Q" ia,x /Q s to be adjusted as a fit 
parameter for each value of q. See figure [T} 

It is well understood that the comparison of saturation ideas to total cross section cr 7 . p 
in DIS data can be made to work with a saturation scale Q s that is independent of impact 
parameter b out to a hard cutoff beyond which it vanishes. The plot with q = 0.1 GeV in 



It is reassuring that the value 6 max = 0.69 fm is close to the root-mean-square charge radius of the proton, 
R p = 0.88 fm, as measured in elastic processes, though perhaps it is a bit surprising to find 6 max < R p . Note 
however that 6 max is sensitive to the overall normalization of cr 7 » p , so it is affected by any uncertainties in 
the overall normalizations of the photon wave-function and of the cross-section. 
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q = 0.1GeV, Qf^ / Q s = 1 




q = 0.2GeV, (2™™ / g s = 1 .2 





q = 0.5 GeV, 




bmax = 18 fm 















*W = 18 fill 




b(fm) 



Figure 1: The plots of a/™Vp show comparisons between predictions derived from (45) for 
various values of 6 max and DIS data, taken from Figure 2 of [211] . In each plot of y/ra 7 * p , 
QT ax /Qs was se t to a constant that led to a good fit with data for 6 max = 0.69 fm. The plot 
at lower right shows the dependence of Q s /Q s 011 b, given by (46), with the values of q and 
QT^/Qs as shown in the other three plots. 
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figure [I] is essentially a confirmation of this point. But, as is evident from the other two values 
of q shown, the data for <t 7 * p is also consistent with significant variation of Q s , at least for 
the SO(3)-invariant functional form ^ supplemented by the hard cutoff in b. Admittedly, 
using q = 0.5 GeV pushes past the envelope where use of (27) and (28) is uniformly reliable 
across the experimentally accessible range of Q. Thus an improved analysis would again be 
desirable. Comparison to diffractive cross-sections would also be helpful. 



5 Conclusions 

The main claims of this paper are contained in equations ^ and ([7]). The first result is the 
statement that there are solutions to the leading order BK equation which respect an 50(3) 
subgroup of the conformal group. The SO (3) subgroup is characterized by a parameter q with 
dimensions of inverse length. SO (3) q - invariant solutions are almost as simple as the more 
widely studied translationally invariant solutions, where the scattering amplitude S depends 
only on the dipole size, not its position in the transverse plane; indeed, translationally 
invariant solutions can be recovered from the q — > limit of S'0(3) g -invariant solutions. After 
mapping the transverse plane stereographically to a sphere, one can understand SO(3) q - 
invariance as uniformity on the sphere. Indeed, the distance function d q appearing in ^ 
and ^ is the the chordal distance between the two images points on the sphere, which is 
obviously 50(3) g -invariant. The form of d q leads directly to the specific dependence of the 
saturation scale Q s on impact parameter shown in ([7]). 

It is noteworthy that in solutions to the BK equation of the form (|6]), the saturation scale 
doesn't change its functional dependence on impact parameter as the rapidity increases. 
Heuristically, the hadronic target doesn't spread out at all, it just gets blacker all over. 

Conformal invariance in QCD is modified by the running of the coupling, and most 
spectacularly by confinement. So we should expect to start to have trouble matching data 
as the characteristic momentum scale approaches Aqcd- Indeed the analysis of section H] 
showed too fast an increase of the cross-section o" 7 * p with r = Q 2 /Q 2 S as r becomes small. 
This is because conformal invariance predicts a fatter tail of energy density and saturation 
scale at large impact parameter than a confining theory like QCD will support. However, 
this is hardly sufficient reason to abandon the conformal approach: when Q™ ax 3> Aqcd 
there may be a substantial region in the transverse plane where the results ^ and (J7]) 
provide a good leading order description of the variation of parton distribution functions 
over the transverse plane. Using conformal invariance as an organizing principle for the 
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early stages of a collision is likely to be especially useful when the early time dynamics itself 
is conformally invariant. This is true in the glasma description [31] (see also the earlier work 
[35] . and [36]), where the relevant dynamics is classical Yang-Mills theory. To the extent 
that conformal symmetry is a good symmetry of both the initial state and the subsequent 
early-time dynamics, the analytic solution of [I] should describe the fluid flow until such 
time as the non-conformal nature of the equation of state becomes important. 

Clearly it would be desirable from a phenomenological standpoint to study deviations 
from SO(3)g invariance in solutions to the BK equation, as I did with A. Yarom for hydrody- 
namics in [SJ. For example, I would like to know whether and how fast the S , 0(3) g -symmetric 
configuration is approached according to the leading-order dynamics 0. More ambitiously, 
one could inquire how NLO corrections cause the inverse size parameter q to evolve with 
rapidity, and how deviations from conformality evolve according to NLO dynamics. Best of 
all would be to develop a quantitative understanding of how deviations from conformality in 
the initial state propagate through to the hydrodynamical regime. 
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